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Low-energy nuclear structure is not sensitive enough to resolve fine details of nucleon- 
nucleon (NN) interaction. Insensitivity of infrared physics to the details of short-range 
strong interaction allows for consistent, free of ultraviolet divergences, formulation of 
local theory at the level of local energy density functional (LEDF) including, on the 
same footing, both particle-hole as well as particle-particle channels. Major difficulty is 
related to parameterization of the nuclear LEDF and its density dependence. It is argued 
that structural simplicity of terminating or isomeric states offers invaluable source of 
informations that can be used for fine-tuning of the NN interaction in general and the 
nuclear LEDF parameters in particular. Practical applications of terminating states at 
the level of LEDF and nuclear shell-model are discussed. 

1. Introduction 

The atomic nuclei are very complex finite many-body systems exhibiting non-trivial 
coupling of single-particle (sp) and collective degrees of freedom. It is well known 
that significant fraction of these many-body correlations can be taken into account 
through symmetry violating intrinsic states within self-consistent mean-field (MF) 
approximation using finite-range Gogny 1 or contact Skyrme 2 interactions. These 
interactions, or the underlying energy density functionals (EDF), arc parametrized 
usually by about ten coupling constants which are adjusted to basic properties of 
nuclear matter and to selected data on finite nuclei. 

The fitting procedure is by no means unambiguous. There is no clear consensus 
of what dataset should be used in adjusting EDF parameters. Coupling between 
single-particle and collective degrees of freedom and correlations beyond MF man- 
ifestly appearing in any finite mesoscopic system implies further that different phi- 
losophy (procedures and datasets) should, in principle, be applied while fitting EDF 
for pure MF applications and for theoretical methods taking explicitly correlations 
beyond MF approximation like random phase approximation, generator-coordinate 
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method, or symmetry-projection techniques to avoid double counting. Yet another 
open problem is related to density dependence of nuclear EDF coupling constants 
and its relation to the NNN interaction. Presently used prescription is purely phe- 
nomenological. It is not only unsatisfactory from theoretical point of view but cause 
numerous troubles in particular in angular momentum projected calculations, see 3 
and refs. quoted therein. Similar ambiguities related to the form and density de- 
pendence apply to pairing channel. 

The number of existing parameterizations of, in particular, Skyrme forces re- 
flects ambiguities listed above. This rather frustrating situation can be healed to a 
certain extent by turning the attention to specific, extremely simple nuclear states 
where pairing, single-particle and collective degrees of freedom decouple to the 
largest possible extent. The classical examples of such states are superdeformed 
states or terminating and high-spin isomeric states, see 4,5 and refs. cited therein. 

Hereafter I will focus on applications of terminating and high-spin isomeric 
states. The aim is to demonstrate how the structural simplicity of these states can 
be used to probe nuclear Skyrme EDF (see Ref. 6 ' 7 for further details), pairing 
correlations as well as effective sdfp shell-model (SM) NN interaction. 

I will start by presenting general arguments speaking in favor of local nuclear 
theory. I will show next a couple of numerical results pertaining to isovector pair- 
ing correlations. In particular, I will consider three versions of density dependent 
delta interaction (DDDI) including volume-active, surface-active and mixed vari- 
ants. I will present examples of calculated shape-gap correlation plots that seem 
to quite firmly point out toward volume-like character of nuclconic isovector pair 
field by excluding both surface-active or mixed pairing scenarios. Next, I will briefly 
overview problems with conventional pair-blocking encountered in the analysis of 
iV=83 high-spin isomeric states in rare-earth nuclei. Finally, I will discuss the energy 
differences between two terminating configurations [/ 7 ™ 2 ]j m „ and [4^/2 /"/t^ An ax m 
Z<N, ^4^45 mass region 

AE = E{[d-} 2 f^] Imax ) £?([# /2 ]/ ma J. (1) 

This extremely simple observable can be used not only to test time-odd spin-fields 
and spin-orbit strength of the nuclear EDF 6,7 but also elucidates much more subtle 
effects directly pertaining to isospin dependence of SM matrix elements. 

2. Infrared nuclear theory 
2.1. Particle-hole channel 

Nuclear structure theory aims to describe low-energy nuclear excitations. Since such 
observables are be definition not sensitive enough to resolve details of the under- 
lying NN (and NNN) interaction one can freely relax any less or more stringent 
relation between effective and ab initio NN (and NNN) potentials and attempt to 
built the effective nuclear theory essentially from scratch using, as the only guid- 
ing principle, the following general statement: low-energy or INFRARED physics 
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should not depend on high-energy or ULTRAVIOLET dynamics. Such separation 
of scales is the underlying principle of effective field theory which aims to describe 
composite objects at low-energies by Largrangians which include ultraviolet dy- 
namics by means of a serie of contact corrections 8 . It implies that the following 
expansion of an arbitrary short-range, rotationally invariant NN interaction: 

vs(q 2 )^vs(0)+v^(0)q 2 +vf\0)q 4 ..., (2) 

in terms of transferred momentum q should converge fast. a In this way complicated 
and in fact unknown two-body NN interaction is mapped by just a few constants: 
vs{0), v^\o), . . . which should be carefully adjusted using representative set of 
nuclear data. 

In the r-space the expansion (2) takes the following form 

v eff {r) w vi ong (r) 
+ ca 2 S a (r) 

+ d ia 4 V 2 5 a (r) + d 2 a 4 VS a (r)V 
+ ... 

+ gi a n+2 V n S a (r) + ... , (3) 

where 5 a (r) denotes an arbitrary model of the Dirac delta function of range a. The 
viong(r) describes the long-range part of the NN potential. Since typical range of 
strong interaction is a ~ 1 fm and is much smaller than nuclear radius R, a <C R, 
the vi ong (r) pertains essentially to Coulomb interaction. The correcting terms are 
proportional to S a (r) and derivatives of S a (r) arranged as in Eq. (3) to assure 
spherical symmetry. Each correcting term introduces its own dimensionless cou- 
pling constant c, d\ , di . . . g\ . . . since the effective range ~ a or ultraviolet cut-off 
momentum <~ 1/a is pulled out explicitly in Eq. (3). As already mentioned these 
coupling constants should be readjusted to a selected set of low-energy data. 

The assumed spherical symmetry defines uniquely (Eq. (3)) the form of the cor- 
recting potential replacing the low-momentum part of the exact potential removed 
by the ultraviolet cut-off procedure. However, since S a (r) can be modeled in es- 
sentially arbitrary way, there is, at least in principle, an infinite set of equivalent 
realizations of the effective interactions or effective theories. In particular, assuming 
Gaussian form factor: 

e -r- 2 /2a 2 

or more precisely modeling i5 a ('') by a sum of attractive and repulsive Gaussians 
of different ranges with space-, spin- and isospin-exchange term and supplementing 
it by a density dependent term and a spin-orbit term (last two terms in Eq. (5)) 

a At this stage we omit, for the sake of simplicity, spin and isospin degrees of freedom. 
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leads to the Gogny force 1 

2 
3=1 

+ t 3 (l + x 3 P a )p^R)6(r 12 ) 

+ iWoicn + <r 2 ) (k' x 5{r 12 )k) ■ (5) 

Since the early work of Brink and Vautherin 9 it is known that particle-hole 
(ph) channel can be described by contact lim a ^o °~a = S(r) Skyrme interaction 2 : 

v(l,2) = t (l + x P a )5(r 12 ) 

+ |t!(l + xiP a ) (k'H{r 12 ) + S(r 12 )k 2 ) 
+ t 2 (l + x 2 P <y )k'S(r 12 )k 
+ ^t 3 (l + x 3 P„) P 2(R)6(r 12 ) 

+ ^VFo(c^ 1 + cr 2 ) (k' x <5(r 12 )fc) , (6) 

where: r 12 = n r 2 ; -R = (ri + r 2 )/2; the momentum operator fc = ^(Vi — V 2 ) 
acts to the right while k' = — ^(Vi — V 2 ) acts to the left hand side. The first 
three terms of the Skyrme interaction correspond to the first three terms in the 
expansion (3). Last two terms in Eq. (6) denote density dependent and spin-orbit 
terms. 

The concept of effective theory based on renormalization of ultraviolet dynamics 
turns up side down the philosophy behind the effective NN forces. It states that due 
to poor resolution of low-energy data the exact form of the NN force is not at all 
needed in practical computations. It can be replaced by local corrections of the form 
given in Eq. (3). The insensitivity of the effective theory to the short-range details 
tells us that we can construct infinitely many theories having the same low-energy 
behavior. All of them (although different and having probably not much in common 
with the true short-range part of the NN interaction) are essentially equivalent in 
the sense that all of them should be capable to reproduce low-energy nuclear data 
with desired accuracy when order-by-ordcr refinement (3) is applied. In this sense 
Gogny and Skyrme interactions appear as two independent (among infinitely many) 
realizations of an effective theory. So far, Skyrme interaction was viewed rather as 
a short-range expansion of finite-range Gaussian force 9 and regarded as a limiting 
case of a seemingly more fundamental finite-range Gogny force. 
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Fig. 1. Mean neutron pair gaps calculated using HFB D1S Gogny method (open dots) in Ref. 10 
in comparison to empirical three-point OES data A (odd — N) (filled dots). 



2.2. Particle-particle channel 

The particle-hole and particle-particle (pp) or pairing interactions can be treated 
independently. One can therefore apply similar expansion to the pp interaction: 

M 2 )~ 9 + 92q 2 +9iq A ■■■ (7) 



u pair \ 



By retaining only the first term in (7) and by modeling it by a Gaussian-type Dirac- 
delta model (4) one obtains finite range Gogny pairing force which appeared to be 
very successful in numerous practical applications. The example of such calculations 
illustrating Hartree-Fock-Bogolyubov (HFB) mean neutron gap calculated using 
D1S Gogny force versus neutron number 10 is shown in Fig. 1. Note, that the 
theoretical pair-gaps follow very closely empirical three-point odd-even staggering 
(OES) data: 



A(JV) 



(-1) 



N 



-[B(N - 1) + B(N + 1) - 2B(N)] 



(8) 



extracted for odd-A^ nuclei which, in accordance with the analysis of Rcfs. n > 12 ; 
represents mean pair gap. In particular no low-mass enhancement of the pair gap 
following conventional textbook A ~ 1/y/A estimate is seen neither in the data 
nor in the calculations. An additional advantage of finite-range pairing model is 
that it automatically removes high-momentum scattering processes. Indeed, due to 
the finite-range r Q (r a ~ 1 fm) the Gogny force discriminates states above E c ~ 
p 2 c /2m r ~ h 2 /mr ~ 40MeV [m r = m/2 is reduced mass] since r Q p c ~ h. 



In spite of its success the use of finite range interaction in the pp channel may be 
viewed as rather unnecessary complication. Indeed, in the nuclear matter pairing 
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Fig. 2. Neutron pair-gaps versus neutron number. Upper left panel shows empirical OES. The 
remaining panels illustrate pair-gaps obtained within Skyrme-HFB model using three different 
versions of the DDDI interaction: surface-active (lower left), mixed (upper right) and volume- 
active (lower right). Taken from Ref. 13 . 



modifies the nucleonic motion essentially only in the closest vicinity of the Fermi 
energy, Ep — A < ^ PF ^^ < Ep + A giving rise to uncertainty in momentum 
space, Sp ~ A/vp. This uncertainty translates to uncertainty in coordinate space 



(hc) A k F 
(mc 2 )A 



3> r Q ~ y- which by far exceeds the typical interaction range r a . Conse- 



quently, the nucleonic Cooper pairs appear as spatially extended objects which can 
hardly feel subtle details of the underlying NN pairing interaction. Nuclear pair- 
ing should be therefore well described within local approximation using plain delta 
interaction or, slightly more general, density dependent delta interaction (DDDI): 



u pair 



(r) 



1 



S(r). 



(9) 



We are again touching a subtle question of resolution and sensitivity of low 
energy data but this time with respect to pairing. In particular case of nuclear 
matter it was demonstrated explicitly by Garrido et al. 14 that one can rather 
easily parametrize interaction (9) in order to reproduce pair-gap A(kp) versus kp 
dependence in nuclear matter (with free particle spectrum) obtained using Gogny 
interaction. In the case of finite nuclei the situation is far more complicated. Figure 2 
illustrates empirical neutron-gap (OES) and the results of large scale spherical 
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Skyrmc-HFB calculations using DDDI interaction in pp channel 13 . Three different 
versions of the calculations are depicted corresponding to surface (p c = p Q ), mixed 
(p c = 2p Q ) and volume (p c — > oo) type pairing scenarios. In all cases a = 1. The 
strength v was adjusted separately for each pairing-mode in order to reproduce 
empirical neutron pair-gap in 120 Sn. 

The overall agreement between neutron OES and the calculated pair-gaps is 
satisfactory for volume-active and mixed pairing scenarios. For these two pairing 
scenarios the agreement is also similar to the one obtained with Gogny-pairing 
shown in Fig. 1. Clearly, the OES itself is not sensitive enough to differentiate 
neither between these two variants of the DDDI interaction nor between the DDDI 
and Gogny pairing. On the other hand the calculations clearly prefer volume or 
mixed DDDI over pure surface-active DDDI. Indeed, the surface-active pairing leads 
to enhanced neutron gaps in light systems which are not at all observed in the data. 

Further inside into a character of nuclear pairing can be gained by looking into 
correlation plots of deformation, which is rather robust observable, versus pairing 
gap. An example of such shape-gap consistency plot is shown in Fig. 3. Two different 
examples of rather well deformed nuclei are depicted in the figure. Left part shows 
/3 2 -vs-A„ plot in 50 Cr. This is a typical case when shape-gap consistency between 
calculations and the data (marked by the area where horizontal and vertical shaded 
areas corresponding to empirical uncertainties in /3 2 and A„ cross each other) can be 
reached essentially irrelevant of the assumed pairing scenario. Hence, this example 
shows no sensitivity at all with respect to the considered pairing variant. The right 
hand side of Fig. 3 shows similar /? 2 -vs-A„ plot but for 46 Ti. In this case shape- 
gap consistency requirement simply excludes both the surface-active and mixed 
pairing scenarios. There are at least two almost obvious questions which need to 
be answered in this context: (i) Can the hierarchy of deformed-to-spherical phase 
transition versus pairing type which so clearly seen in Fig. 3 be understood in a 
simple manner? (ii) Is the case exclusiveness in 46 Ti accidental or it can be traced 
down systematically throughout the periodic table? Studies along these lines are 
under way. 

2.3. Toward consistent superfluid local density approximation 

The disadvantage of using DDDI type (local) pairing interaction as compared to 
finite-range Gogny interaction is that the former requires explicit cutoff parameters 
to avoid divergences. In particular, all Skyrme-HFB calculations presented above 
were done within limited phase-space for pairing. Although practically sufficient 
in most applications, such a brute-force method is unsatisfactory from theoreti- 
cal point of view. In this respect, local pairing can be considered as considerably 
disadvantageous as compared to finite-range pairing. It appears however, that the 
divergence can be rather easily identified and subsequently regularized leading to 
consistent cutoff-free superfluid local density approximation (SLDA) 15 ' 16 . 

The appropriate scheme for regularization of ultraviolet divergence in anomalous 
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Fig. 3. Shape-gap correlation plot in 50 Cr (left) and 46 Ti (right). Three different curves represent 
Skyrme-HFB calculations assuming surface-active (open dots), mixed (triangles) and volume- 
active (black dots) pairing variants. 



density matrix: 



v(r u r 2 ) =^<(riK(r 2 ) ~ _ , (10) 



at the level of the local density approximation (LDA) i.e. including properly dom- 
inant particle- hole channel was recently proposed by Bulgac and Yu 16 . The idea 
is to renormalize divergent terms by introducing cutoff (E c = ) dependent 

counter-terms. Such an approach leads to standard local HFB formalism with cut- 
off parameters but with a gap equation dependent on the effective running coupling 
constant: 

E c 

v c (r)=Y / v*(r)u i (r), (11) 

Ei>0 

A(r) = -g eff (r)u c (r), (12) 
1 1 m{r)k c {r)\ i k F {r) ^ k c {r) + k F {r) j ^ 



9eff(r) g[p(r)} 27r 2 fi 2 1 2k c (r) k c (r) - k F (r) 

Introducing a running coupling constant implies that the cutoff dependence is only 
formal and disappears for sufficiently large E c 16 . The cutoff- free superfluid LDA 
approach is now in phase of extensive tests 16 ' 17 . 

3. Pairing in high-spin isomeric states 

The atomic ground states are strongly correlated. Simple mean-ficld-plus-pairing 
model, although quite successful in reproducing nuclear masses, see Ref. 18 and refs. 
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cited therein, cannot incorporate all important correlations and polarization effects 
within single symmetry- broken Slater determinant. Hence, in order to test/tune less 
ambiguously basic theoretical ingredients of our models is seems natural to change 
the attitude and look not into terribly complicated strongly correlated states but 
turn the attention toward as pure as possible physical situations. 

The high-spin isomeric states (HSI) open new and so far not fully explored possi- 
bilities to study, in particular, pair correlations, blocking phenomena and superfluid- 
to-normal phase transition 19 > 20 . They are structurally extremely simple hence both 
configuration and shape can be kept rather well under control in the calculations. 
In turn, shape and pairing polarization due to blocking can be studied in detail. 
As an example let us recall systematic calculations by Xu et al. 21 who, using 
HSI, re-examined traditional average gap method 22 used to determine monopole 
pairing strength Gmn for Lipkin-Nogami calculations. In particular, it was found 
that inclusion of polarization effects requires ~10% larger pairing strength Gmn 
as compared to the average-gap value 21 in order to reproduce experimental data. 

Recently the HSI have been systematically observed in N=83 nuclei with 
60<Z<67 23 . In odd-A nuclei the HSI have J w =49/2 + and correspond to a seniority- 
five, stretched shell-model configuration [C : ^(f7/2hg/2ii3/2)<8> 7i"h 11 / 2 2 ]49/2 + - In 
odd-odd nuclei the observed HSI have J 7r =27 + and are assigned to [C DO : 
^(f7/2hg/ 2 ii3/2)<8> 7i'(d^, 1 2 h 11 /2 2 )]27 + configuration. These two configurations differ 
by a single proton hole in Nilsson [402] 5/2 orbital originating from spherical d$/2 
sub-shell. This unique data set enables to study for the first time OES both at 
the ground states (GS) as well as at high-spins using conventional technique based 
on binding energy indicators 24 . The most striking feature of this data set is an 
almost constant excitation energy of the HSI what implies that OES at the ground 
states (GS) follows very closely the HSI value Ags{Z) ~ Ahsi{Z). Conventional 
interpretation of this empirical result in terms of pairing-gap suggests the lack of 
blocking phenomenon. Note, that this conclusion is independent on the type of 
pairing indicator. Hence, in the following, we will use standard three-point filter 
A(Z) (see 11,12 for detailed discussion of its physical content) of Eq. (8). 

The importance of pair-correlations at the HSI can be elucidated by considering 
extreme single-particle model. An example of such calculations is shown in Fig. 4. 
The figure compares experimental OES of Eq. (8) to theoretical OES computed 
using self-consistent SHF model with SLy4 25 and SkO 26 parameterizations. The 
theoretical values show clear systematic rise not at all observed in the data. This 
theoretical trend is generic and reflects the fact that, within the extreme sp sce- 
nario, the indicator (8) measures simply a distance between proton Fermi energy 
ep and proton hole state e[ 40 2]5/2 as indicated in the right hand side of the figure. 
The increase of A(Z) with Z is therefore generic to any sp model. It can be natu- 
rally stopped (or slowed down) when ph excitations are replaced by quasi-particle 
excitations i.e. in the presence of relatively strong pair-correlations. 
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Fig. 4. Experimental (dots) and theoretical (squares) OES of Eq. (8). Calculations for fixed con- 
figurations C and C 00 were carried out with single-particle SHF-SLy4 and SHF-SkO models. 
Right hand side shows spherical single-particle spectra of all models discussed in the context of 
HSI normalized to Tin/2- This part of the figure illustrates schematically a dominant within the 
sp scenario and increasing with Z contribution to OES, A(Z) ~ ep — e[ 4 Q 2 ]5/2, originating from 
energy difference between "fixed in energy" proton hole [402] 5/2 and proton Fermi energy which 
is "moving up" in energy with increasing Z. 
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Fig. 5. Experimental (filled circles) and theoretical (DfPM: squares; DPES: triangles) excitation 
energies of HSI in the N=83 isotones. Open triangles show the DPES calculations corrected by 
the pn residual interaction 8 pn extracted from nuclear binding energies using the 9-point indicator 
of Ref. 24 . 



To study pairing properties of GS and HSI states we employed two different 
approaches: the deformed independent particle model (DIPM) 27 and the diabatic 
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Fig. 6. A(Z) of Eq. (8) calculated using DPES model. Filled (open) symbols denote GS (HSI) 
values. The effect of S pn on GS DPES values is indicated (gray triangles). Experimental values of 
A(Z) are marked by dots. 

potential energy surface model (DPES) 20 . For the 7V=83 isotones, both DIPM and 
DPES methods yield fairly consistent results. Namely, both methods predict: (i) 
weakly deformed ground states; (ii) well-deformed, oblate (/3 2 ~ —0.2) HSI states 
(in) C and C 00 yrast HSI configurations. The only exception is 150 Ho where DPES 
predicts the C QO configuration to lie ^150 kcV above the HSI configuration involv- 
ing four aligned hn/ 2 protons. The calculated excitation energies of HSI, AEhsi, 
slightly depend on the model used reflecting mostly the differences in sp spectra 
of the underlying mean-potentials, see Fig. 4 and Fig. 5. Note, that agreement 
between DPES prediction and experimental AEhsi is excellent particularly after 
correcting theoretical results for residual neutron-proton interaction present in the 
GS of odd-odd nuclei. Further details concerning both the models and results can 
be found in 23 . 

The OES of Eq. (8) calculated using the DPES model is shown in Fig. 6. The 
values calculated at the HSI follow experimental data very closely meaning that our 
theoretical tools and/or interactions can be very reliably tuned out in such simple 
situations. On the other hand, the theory faces serious troubles in reproducing OES 
in the GS. The complexity of GS causes that contributions to OES due to pairing, 
sp-proton energy splitting, residual proton- neutron (pn) interaction and/or nuclear 
magnetism (time-odd effects) which must all be taken into account simultaneously 
in order to reproduce experimental data are difficult to control in a satisfactory way 
in strongly correlated states. In particular, our calculations clearly indicate that 
contribution due to blocking which is mainly responsible for shift of theoretical GS 
and HSI curves in Fig. 6 is far too strong and requires revisiting. 
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4. Probing isospin dependence of SM matrix elements at band 
termination. 

Recently, it was shown 6 that a set of terminating states in the mass region 
45 may provide a unique tool to constrain the Skyrme energy density functional 
(SEDF). The key idea is to look into energy differences (1) between two types 
of terminating configurations [fy2\i max and [d^, 2 f"i ]i max i- c - mto tnc quantity 
which, as the long experience in the studies of terminating states within the Nils- 
son model have clearly shown, see Ref. 4,5 and refs. quoted therein, belong to the 
purest single-particle observables available in nuclear structure. In particular, it was 
demonstrated 6 that by constraining the SEDF to the empirical spin-isospin Lan- 
dau parameters and by slightly reducing the spin-orbit strength, good agreement 
with the data could be obtained. This result, based on high-spin data for terminat- 
ing states, is consistent with conclusions of previous works 28 ' 29 based on different 
theoretical methodology and experimental input (such as giant resonances, beta 
decays, and moments of inertia). The validation of the assumption 6 regarding the 
single-particle character of the maximally-aligned states and the robustness of ter- 
minating states in determining properties of the nuclear EDF was further supported 
by the recent comparative study between the fully correlated shell model (SM) and 
the Skyrme Hartree-Fock (SHF) 30 , showing essentially a one to one correspondence 
between the models at least for N>Z nuclei, see Fig. 7. 

In the N=Z nuclei the comparison between the models can be done only after 
evaluation of the correlation energy SEt due to the spontaneous breaking of isobaric 
symmetry in the SHF. The correction due to isospin symmetry restoration shifts 
the T=0 state down in energy in the laboratory system, as illustrated in the inset of 
Fig. 7. After (phcnomenological) symmetry restoration predictions from both SM 
and SHF models become very similar in all considered nuclei. However, while the 
agreement in N>Z is satisfactory the N—Z results deviate strongly from the data. 
Definitely too strongly as compared, in particular, to overall excellent performance 
of the fp shell SM 31 . 

One possible origin of this deviation may have its source in the assumed trun- 
cation to lp-lh cross-shell excitations. This configuration-space restriction is ex- 
pected to impact the isoscalar channel associated with the sd-^fp pair scattering 
i.e. isoscalar proton-neutron pairing (pn) mode which is also absent in conventional 
mean-field description. A possibility of an onset of the isoscalar pn pairing nearby 
band-termination was already discussed within the mean- field formalism in Ref. 32 . 

An alternative source of deviation may be traced back to incorrect isospin de- 
pendence of SM matrix elements. The quantity (1) involves particle-hole excitation. 
In such a case, according to Bansal and French 33 , and Zamick 34 (BFZ) ... it is a 
grave error to assume that ph interaction is independent on isotopic spin.... Under 
certain assumptions listed in Refs. 33 > 34 the isotopic dependence of ph interaction 
can be approximated by a simple monopole interaction <~ bt\ - t<i yielding additional 
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Fig. 7. Difference Ai?th — Ai? cxp between experimental and theoretical values of AE (1) in A~44 
mass region. Dots denote the SM results. Circles denote the SHF results based on the modified 
SkO parameterization (see text). The SHF calculations for the \d~^\S™t^\l max intruders in N=Z 
nuclei yield two nearly degenerate states associated with proton (7r p h) and neutron (i-'ph) cross- 
shell excitations. As shown in the inset, the physical T=0 state in the laboratory frame is shifted 
down in energy by 5Et (isospin correlation energy). Squares denote the SHF results for N=Z 
nuclei with the isospin correction added. The SHF results were shifted by 480 kcV in order to 
facilitate the comparison with SM. Taken from Ref. 30 . 



contribution to AE of Eq. (1) of the form: 

AE T = h(T(T + 1) - T p (T p + 1) - T h (T h + 1)), (14) 

where Th = 1/2 and T p = T ±1/2 are hole and particle contributions to the 
isotopic spin T. In the case of terminating states in N=Z nuclei T p = 1/2 while 
for terminating states in N>Z we have T p = T — 1/2. Hence, the anticipated 
contributions to the AE due to the BFZ mechanism are: 

!--b in N = Z 

1 / 1\ as) 
- 2 b[T--) in N*Z 

If, for some reason, the strength of the monopole interaction b is overestimated in 
the SM by 56=700 keV (in this mass region it amounts to 5b/b^20%) the SM values 
of AE should be corrected by: 525 keV in N=Z nuclei; OkcV in T = 1/2 43 Sc, 45 Ti, 
and 47 V nuclei; -175 keV in T = 1 42 Ca, 44 Sc and 46 Ti nuclei; -350 keV in T = 3/2 
nucleus 45 Sc; — 525 keV in T = 2 nucleus 44 Ca. The corrected values are depicted 
in Fig. 8. Note that they match empirical the data almost perfectly. 

The suggestion that the discrepancy between the SM description of N > Z and 
N = Z nuclei may be caused (fully or partly) by incorrect isospin dependence of 
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Fig. 8. The differences ABjm — AEgxp calculated using SM model (black dots). Grey dots 
indicate SM values corrected by the anticipated contributions due to BFZ effect of Eq. (15), 
assuming <56=700 keV decrease in the b strength. 



the SM matrix elements 35 > 36 is extremely intriguing and the anticipated accuracy 
of AEsm — AEex p deduced solely on the basis of the BFZ formula is indeed very 
appealing, see Fig. 8. To verify the concept we performed two set of SM calculations 
using modified SM interactions marked lateron as SM^ and SM^ 2 ), respectively. 

In the SM« run ALL diagonal particle-hole matrix elements of original SM 
interaction 35,36 were renormalized according to the following scheme: 



V phph (JT = 0) 
V phph (JT=l) 



V phph {JT = 0) + 3d 
V phph {JT = 1) - d 



(16) 



with d = 0.175keV. The use of transformation (16) can be partly justified within sin- 
gle j-shell shell-model phenomenological mass formula where it appears to leave in- 
variant all contributions to the binding energy except the symmetry energy strength 



bsym which changes according to the following rule: b. 



sym 



J sym 



Ad, see Ref. 



37 



It should be underlined, however, that b and b sym although having certain com- 
mon features are different quantities. Such an attempt to modify b was not fully 
successful and satisfactory. On one side we obtained substantial improvement of 
agreement to the data for low-lying particle-hole excitations as compared to the 
original interaction as shown in Fig. 9. On the other side the changes appeared 
not neutral with respect to other observables spoiling, in particular, the previously 
obtained nice agreement between SM and experimental binding energies. 

In the SM< 2 ) run we have decided to apply renormalization scheme (16) to 
CROSS-SHELL matrix elements only. Preliminary calculations shows that such a 
simple procedure leads to an almost perfect agreement between theory and exper- 
iment for both terminating as well as low-lying particle-hole excitations. Further 



I 



I 
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Fig. 9. Difference between the SM and experimental excitation energies for selected low-spin lp- 
lh intruder states. Calculations done using original SM interaction (bullets) are compared to the 
results obtained for modified SMW interaction. Note dramatic improvement obtained for modified 
SM. 

details concerning our systematic refinements of SM interaction inspired by very 
precise high-spin data including both sets of modifications discussed above will be 
published elsewhere 38 . 

5. Summary and outlook 

It is shown, that local effective theory suitable for low-energy nuclear structure 
calculations can be consistently formulated as free form ultraviolet divergences su- 
pcrfluid local energy density functional approach. The underlying assumption con- 
cerning independence of infrared physics on ultraviolet dynamics implies that such 
a strictly local theory should in principle be equivalent to finite-range effective the- 
ory concerning accuracy of its predictions. This conclusion makes local realization 
of effective theory extremely appealing due to its numerical simplicity. 

There are at least two major problems related to the nuclear LEDF including: 
(i) Density dependence of LEDF parameters or coupling constants and its relation 
to effective three body NNN interaction. This problem pertains both to local as 
well as to finite-range nuclear EDF formalism in the same way. (ii) Datasets used 
to parametrize are definitely incomplete leading to a multitude of local (Skyrme) 
effective interactions. In this respect finite-range realizations of effective theories are 
less prone to the details of datasets as our practice with the Gogny interaction shows 
where essentially only two parameterizations are used in numerical applications 39 . 
This is related with smaller number of terms that are used in the expansion (3). 

The relations between the LEDF parameters and fine-tuning of the LEDF pa- 
rameters can be reliably done in simple physical situations. Long standing experi- 
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ence tells that terminating or isomeric states belong to the purest known examples 
of almost unperturbed single-particle motion 4 . In this work we summarize our re- 
cent efforts related to terminating and isomeric states showing, in particular, that: 

• The structural simplicity of these states can be used to unify time-odd 
spin-fields and tune up spin-orbit strength of the LEDF, see 6 ' 30 for further 
details. 

• They can be used for identification, evaluation and subsequent restoration 
of broken symmetries inherently obscuring the SHF treatment, see Rcf. 30 - ' 
for further details. 

• Multi quasi-particle configurations in rare-earth nuclei can be used to test 
pairing interaction and blocking phenomena 23 . 

• Finally, it is shown that these states appear to be extremely useful in iden- 
tifying and correcting isobaric dependence of cross-shell matrix elements of 
sdpf SM interaction 30 ' 38 . 
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